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Abstract. We present simulations of initially stable isothermal clouds exposed to ionizing radi- 
>D ation from a discrete external source, and identify the conditions that lead to radiatively driven 

" ' implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN 
' (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas 
f^ , et al. (2009). We find that the incident ionizing flux is the critical parameter determining the 
j.^ ' evolution: high fluxes simply disperse the cloud, whereas low fluxes trigger star formation. We 

flnd a clear connection between the intensity of the incident flux and the parameters of star 
yS^ , formation. 

Oh- 

' ' Keywords, hydrodynamics, methods: numerical, stars: formation, (ISM:) HII regions 



o 



►^ , 1. Introduction 

C . When an expanding Hii region overruns a pre-existing cloud, it compresses it by driving 

^^ ' an ionization front and a shock wave into it (Sandford et al. 1982; Bertoldi 1989; Lefloch & 
p^ Lazareff 1994). The inner parts may become gravitationally unstable and collapse to form 

new stars. This mechanism is known as Radiation Driven Implosion (RDI). Observations 
(Lefloch & LazarefF 1995; Lefloch et al. 1997; Sugitani et al. 1999, 2000; Ikeda et al. 
^^ ' 2008; Morgan et al. 2008; Chahuan et al. 2009) strongly support a connection between 
T— I the RDI mechanism and the formation of Young Stellar Objects (YSO). Simulations of 

s/ the interaction of the ultraviolet ionizing radiation with self-gravitating clouds have been 

• ^H . presented by various authors (Kessel-Deynet & Burkert 2003; Gritschneder et al. 2009; 
^\ ' Miao et al. 2009). However, no model can explain where star formation takes place (in 
the core or at the periphery) or when (during the maximum compression phase or earlier 
- Deharverg et al. 2005). 

In this paper we perform a set of 75 simulations of clouds exposed to ionizing radiation. 
The aim of this work is to answer questions of whether the incident ionizing flux is able 
to trigger the formation of new stars or not, and how the process and the properties of 
this star formation are connected with the intensity of the incident flux. In Section 2 we 
give a brief description of the numerical treatment and the initial conditions we use. In 
Section 4 we discuss the results of our simulations. We summarize in Section 5. 
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2. Numerical Treatment and Initial Conditions 

We use the Smoothed Particle Hydrodynamics (SPH) code SERE]N(j, fully described 
in Hubber et al. (2010), with an ionization routine (Bisbas et al. 2009) based on the 
HEALPijy sphere tesselation code (Gorski et al. 2005). The temperature of the neu- 
tral gas at density p is T^{p) = T^^^ \} + (p/Pcrit)^~ /> where T^^^ = 10 K, p^mt = 

lO^^'^gcm"^ and 7 — 5/3 is the ratio of specific heats. The temperature of the ionized 
gas is taken to be Ti = 10** K, except in the transition zone between the two extremes, 
where it changes smoothly from Ti to T„ (see Bisbas et al. 2009). We include sink particles 
(Bate et al. 1995) with radii R^iuji — 2.5 AU created if p > p^-^^^^ — 10^^^ gcm^"^. 

Our clouds are stable Bonnor-Ebert spheres (Bonnor 1956; Ebert 1957; heareafter 
'BES') with dimensionless cut-off radii ^3 = 4,5,6 and with masses M = 2,5,1OM0. 
The particle resolution we use is 5 x 10^ SPH particles per solar mass (c.f. Hubber et 
al. 2006). We place the BESs at distance D = lOi?, where R is the radius of the cloud 
(in pc), in order to keep constant the divergence of the incident flux and as parallel as 
possible. We use a single source emitting Lyman-a photons. We run simulations with a 
wide range of emission rates Af^ ^ — 10^ s^^, where x = 48, 48.5, ...52. 



3. Results 

In Fig[T^ we present a semi-logarithmic diagram where we correlate the intensity of 
the incident ionizing flux with the initial mass of each BES. The lines define subsets 
of parameter space where models either show star formation (left) or not (right), with 
accuracy 0.25 dex. It can been seen that as the mass of the BES decreases (and as a 
result ^3 increases) the clouds appear to survive longer in higher fluxes. This is because 
for a given ^3 , the density Pc at the centre of each BES increases with decreasing M . 

We also find that the Str0mgrcn radius at the end of the i?-type expansion determines 
whether stars are formed or not; if the ionization front has not overrun the central core of 
the BES, then the incident flux will trigger star formation during the Z)-type expansion 
of the Hii region. 

The time, ^sinkj between the beginning of the D-type expansion and the first sink 
creation (beginning of star formation) increases with decreasing ionizing fiux. This finding 
is in agreement with simulations of the RDI performed by Gritschneder et al. (2009). 
Figure [TJd is a logarithmic diagram where we plot the values of igmK versus the incident 
flux $j^ . Remarkably, the age of the cloud when star formation occurs does not depend 
on its properties (i.e. the initial mass, M, or ^3). Results from our simulations can be 
described with a power law of the form ^31^^ = 80 x <I)~0-3 {t^^^^ in Myr, <I>j-, in cm^^s^^). 

Figure [2] shows column density plots of a BES with M = lOM© and with f 3 = 6 
at igjNK for different fiuxes. A common feature in all our simulations is that stars form 
close to the symmetry axis joining the centre of the cloud to the exciting star. This is in 
an agreement with observations by Sugitani et al. (1999). The distance dt between the 
first sink particle and the ionization front is a function of the ionizing fiux and the BES 
parameters (see FiglTJ; where we plot dt/2R for all BESs with ^3 =6). We find that for 
low fluxes stars tend to form in the innermost part of the filamentary structure, whereas 
for high fiuxes stars tend to form at the periphery of the cloud. Similar results are found 
also with f 3 = 4 and ^3 = 5. 

Figure [5] shows that the lateral compression, Wd, of the BESs at the beginning of star 

f http://www.astro.group.shef.ac.uk/seren 
X http://healpix.jpl.nasa.gov 
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formation is connected to the intensity of the incident flux. We see that for low fluxes, Wd 
is quite high and the cloud has a U-shape structure, whereas for high fluxes w^ is small 
and the cloud has a V-shape stucture. In FiglTJi we plot Wd/2R for all BBSs with ^3=6 
and we find that stars tend to form during maximum compression once the incident flux 
is increased. Similar results are found also for the rest of the clumps. 



4. Conclusions 

We present simulations of Radiation Driven Implosion in stable clouds represented 
by Bonnor-Bbert spheres. We performed 75 simulations with clouds of different masses, 
different dimensionless radii, and with a wide range of incident fluxes. In general we find a 
connection between the incident ionizing flux and the resultant star formation efficiency. 

We introduce a semi-logarithmic diagram (flux-mass diagram) where we correlate the 
intensity of the incident flux and the initial mass of each BBS, and we deflne zones of Star 
Formation and no-Star Formation. We flnd that if the Str0mgren radius at the end of 
the i?-type expansion has not overrun the central core of the BBS, the ionizing radiation 
will trigger star formation. The age of the cloud when star formation occurs increases 
with decreasing flux, and it does not depend on the properties of the BBS. A power-law 
of the form t^-^^-^ = 80 x $~°-3 flts very well with the results of our models. Finally, as 
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Figure 1. (a) The flux-mass semi- logarithmic diagram where we define areas where stars are 
formed (SF) and areas where stars are not formed (no-SF) depending on the dimensionless 
radius ^g of a BES. (b) Logarithmic diagram of the incident flux versus the age of the cloud, 
^siNKi when star formation occurs. The power law we propose (sohd line) flts very weU with our 
simulations (tsi^K is in Myr and $□ is in cm~^ s~^). (c) Star formation occurs at the periphery 
with increasing flux, (d) Star formation occurs during maximum compression with increasing 
flux. 
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Figure 2. Column density plots of a BES with M = 10 Mq and ^3 = 6 at tgi^K when it is 
exposed to three different intensities of flux (flux increases from left to right). The white dots 
are sink particles. In the left plot we draw the values of di and w^. 

the incident flux increases, stars tend to form closer to the periphery of the cloud and 
during its maximum compression phase. 
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